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We regularise the 3d model by discretising the Euclidean time 
and representing the spatial part on a fuzzy sphere. The latter involves 
a truncated expansion of the field in spherical harmonics. This yields a 
numerically tractable formulation, which constitutes an unconventional al- 
ternative to the lattice. In contrast to the 2d version, the radius R plays 
an independent role. We explore the phase diagram in terms of R and the 
cutoff, as well as the parameters and A. Thus we identify the phases 
of disorder, uniform order and non-uniform order. We compare the result 
to the phase diagrams of the 3d model on a non-commutative torus, and 
of the 2d model on a fuzzy sphere. Our data at strong coupling repro- 
duce accurately the behaviour of a matrix chain, which corresponds to the 
c = 1-model in string theory. This observation enables a conjecture about 
the thermodynamic limit. 
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1 Introduction 

A variety of approaches to the regularisation of quantum field theory exist. 
Dimensional regularisation [1] is most popular in the framework of per- 
turbation theory. In order to overcome the limitations of a perturbative 
expansion, however, a regularisation should restrict the formulation to a 
finite set of degrees of freedom. If the Euclidean action is real and bounded 
from below, a model can then be treated numerically as a statistical sys- 
tem. This method provides in many cases the only access to observables 
beyond perturbation theory or semi-classical approximations. 

From a general perspective, field theoretic models start from some al- 
gebra A for functions on a manifold Ai, and a differential operator V with 
its Hilbert space H. The standard approach for non-perturbative studies 
discretises the manifold to a lattice. Then the degrees of freedom to work 
with are usually the field variables on the lattice sites or links (see e.g. Ref. 



2 



[2]). As an alternative, also Monte Carlo simulations employing the fields 
at discrete momenta have been suggested [5], though much less explored^ 
In both cases A4 is reduced to a finite lattice. 

Generally the goal is to approximate a triple {A, 'H,T>) [4j . This might 
be achieved in quite abstract ways, but in practice a physical picture for 
the regularised system is a useful guide-line. In the lattice formulation 
one approximates the entire triple. The algebra is approximated by a 
commutative algebra on a lattice of points which discretise M., while the 
differential operator is obtained by a finite difference approximation to P, 
and the Hilbert space is adapted to this operator. 

Here we are concerned with an alternative scheme, which is endowed 
with a physical picture on the regularised level as well. Instead of the dis- 
crete eigenvalues of the space-time or momentum coordinates, we now deal 
with angular momentum coordinates. To this end the fields are wrapped 
on a sphere and expanded in spherical harmonics. A related idea occurred 
already in an early construction of a non-commutative space, which added 
an extra dimension and preserved 5d Lorentz symmetry [S]. This method 
benefits from the natural discretisation of angular momentum space in 
quantum physics, but a cutoff still has to be imposed. In the interpreta- 
tion of the angular momenta as spherical coordinates, the cutoff renders 
the sphere fuzzy. These coordinates are embedded in matrices, which are 
Hermitian in the case of the neutral scalar field to be considered here. 

The concept of a fuzzy sphere regularisation has been established in 
Refs. |6]. The Laplace-Beltrami operator T)^ is given in terms of angular 
momentum operators, which are expressed by A^-dimensional irreducible 
representations of SU{2). The algebra A is expanded in the polarisation 
tensors, which are matrix analogues of the spherical harmonics [7j. In 
contrast to the lattice, this regularisation does not explicitly break the 
space symmetries. Its analytic properties have been studied extensively in 
recent years, but the applicability in numerical simulations is less explored. 
The questions are if such simulations are feasible and to what kind of limits 
the measured observables can be extrapolated. So far the 2d A0*^ model has 
been investigated in this respect [SjH Here we extend this study to three 
dimensions, where the spatial plane is mapped onto a fuzzy sphere, and 
the Euclidean time is lattice discretised. This extension entails qualitative 
differences, which are essential in view of the prospects of proceeding to four 
dimensions. In particular the radius R of the sphere plays an independent 
role (it cannot be absorbed by simple rescaling). The recovery of a fiat 

"'^Of course, the standard Hybrid Monte Carlo algorithm for dynamical fermions in- 
cludes a Langevin ingredient in momentum space, but in that case the basic regularisa- 
tion is nevertheless a space-time lattice. 

^Further numerical studies on the fuzzy sphere address U{1) gauge theory j^]. 
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space without truncation requires the hmits N, R ^ oo. 

Section 2 presents the fuzzy sphere formulation of the model, along 
with suitable order parameters. The identification of the phase transitions 
in the (A,m^)-plane is described in Section 3. Section 4 discusses the scal- 
ing of the phase transition lines in terms of N and R. Section 5 compares 
our results to the phase diagram of the corresponding 2d model on a fuzzy 
sphere, and to the model on a 3d non-commutative torus. We also demon- 
strate that our data at strong coupling agree with the behaviour of a matrix 
chain model, which attracted interest in string theory. This observation al- 
lows for a conjecture about the large N extrapolation, as we point out in 
Section 6. Our results are summarised in Section 7, and technicalities of 
the simulation are added in an Appendix. 

A synopsis of this work has been anticipated in a proceeding contribu- 
tion [lOj, and details are presented in a Ph.D. thesis jl 1 j . 



2 The fuzzy sphere formulation of the 3d 
model 

2.1 Regularisation 

In this subsection we specify the regularisation that we used in our sim- 
ulations. The theory to be regularised is the model in 3 dimensions, 
where we assume periodic boundary conditions in the Euclidean time t , 
and the space is taken as a sphere in IR^. Thus the action reads 
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T is the temporal periodicity and R is the radius of the sphere, which we 
are going to denote by 5*^ . (f){t, x) G IR is a scalar field and = Yl^=i ^ 
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Li being the angular momentum components. 

For the regularisation in time we introduce At equidistant sites and 
replace dt by the standard lattice derivative. We are going to use lattice 
units, i.e. we set A'j = T. Thus a configuration is given by a set , 
t = l,...,Ar,. 

Our regularisation of the sphere S\ is less standard, but it also relies on 
a concept established in the literature The coordinates Xi are replaced 
by operators Xj, which still obey the constraint 

3 

Y.Xl = R^-\. (2.2) 
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A truncation to a maximal angular momentum fmax nieans that the oper- 
ators Xi take the form of N x N matrices with = fmax + 1, 

Xi = -j^^= U , {Xi G Matjv) . (2.3) 
V A"' — 1 

The Li are generators in an A^-dimensional irreducible representation of 
SU{2). These coordinate operators do not commute, 

[X^.X,]=i-J^=e,,uXu. (2.4) 

Thus they cannot describe sharp points; the sphere becomes fuzzy. 

A scalar field, which can be expressed as a power series in the coordi- 
nates, now turns into an expansion in the operators Xi (at some fixed time 
site t). Thus this formulation represents the field by A^ x A^ matrices 
In particular for the neutral scalar field these matrices are Hermitian. Its 
spatial derivatives are given as commutators, di4>t{x) — > i [Lj, $f] . 

In summary, the recipe for the regularisation from a sharp to a fuzzy 
sphere involves the replacements 

XiE S\ XiE Mat AT 
04 (x) G C°°(5^) ^ G Mat TV (Hermitian) 

3 
i=l 

[ dnM^,^) ^ T7Tr($0, (2.5) 



where the last relation preserves the normalisation {dQ = smd dd dip). 



We implement these transitions at each discrete time site t. This leads 
to field configurations given by $ = $2, • • • , ^Nt}i and to the action 

Sm = ^Tr[^{i($,,,-$,)^ + ^$,£^$, 

In this regularised form the functional integral reduces to an integration 
over the independent elements of the Hermitian matrices. This is in fact 
tractable in Monte Carlo simulations; note also that the action (12.61) is real. 

A qualitative difference from the 2d model on a fuzzy sphere (which 
corresponds to our model on a single time site) is that the radius R plays 
an independent role; it cannot be absorbed in the coupling constants. 
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A virtue of this approach — compared to the usual space discreti- 
sation — is that continuous spatial rotational symmetry persists on the 
regularised level. The fuzzy sphere is rotated by the adjoint action of an 
element U G SU{2) in the iV-dimensional irreducible representation, 

X U^XU = nX , <l>t ^ U^^tU , (2.7) 

where U can be written in the form U = exp(ia;L) , and TZ G 5*0(3). A 
global rotation in all time sites leaves the action (12. 6p invariant. 

This virtue may prove particularly powerful in cases where continuous 
rotational and translational symmetry (which we obtain in the large R 
limit) play a central role, such as supersymmetric modelsjf] 

The question how profitable that symmetry ultimately is has to be in- 
vestigated based on non-perturbative results. A prerequisite is a controlled 
large N limit, and testing this property is a goal of the current work. It 
should also illuminate the status of possible pitfalls. In particular, the non- 
commutativity of the operators Xi implies a non-locality of the interaction 
in the regularised model. We are going to see that this property can indeed 
affect the thermodynamic limit Q A further goal is to elaborate links of the 
observed universality class to other models of interest. 



2.2 Observables 

Now we introduce the observables to be measured numerically. For this 
purpose we first perform a field decomposition, which is compatible with 
the rotational symmetry. 

The original field (j) can be decomposed in the basis of spherical har- 
monics Yim on S"^ , 

oo e 

0(t, ^,^) = Y^Y, cem{t)Yem{v, i3) . (2.8) 

1=0 m=-l 

In full analogy, the regularised space MatTv has a basis consisting of the 
polarisation tensors Yim, see e.g. Ref. For i = 0, . . . , N—1, m = —£, . . .£ 

■^Literature on supersymmetric systems on a fuzzy sphere exists regarding the theo- 
retical basis [12j and first simulations 13J, though there are many outstanding issues in 
that field. Another important point in this context is that — in addition to the space 
symmetries — also chiral symmetry is intact on the fuzzy sphere, without a fermion 
doubling problem 14J. 

^The impact of a non-local regularisation on the continuum limit is intensively dis- 
cussed in the lattice community (see e.g. Refs. [T5|) in particular in the light of recent 
large-scale QCD simulations with "rooted staggered fermions" . 
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these are N"^ matrices with the characteristic properties 

47r ^ ^ ^ 

^e\-m) = (-l)"^<?m , 

jC%^ = i{i+l)Y,m. (2.9) 
Their construction is reviewed in Ref. [llj. The leading examples are 

^00 = f— ^Af 



with (La),^ = ^{N + 1~ 2%) 



We perform this decomposition in each time site, 

N-l I 



^t = Y.Y. ' (2-10) 



1=0 m=-e 

so that $t is fixed by the A^^ coefficients 

An 



Q™(t) = ^Tr(Fi$,) . (2.11) 
In this work we consider the time averaged terms 

^■=^^'^t, Q„ := -^^Q„(t) , (2.12) 



^ t ^ t 



which are related as 

^47r„ 47r 



Coo = ^Tr($) , ci^ = — Tr(Fil„$j , etc. (2.13) 
We further introduce the quantities 



oo 



l^ll' ■■= E^' = 77Tr($^)- (2-14) 

£=0 
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phase 




disordered 




uniform ordered 




non-uniform ordered 





Table 1: The respective magnitudes of the order parameters {ipo) and {(pi) 
for the three phases that we observed in our results of Section 3 and 4- 
Here the terms in definition ^2.14 ) are approximately related as (||5>||^) k, 
(fo) + (v^i)j ^^^^ higher contributions (fj^i) are small. In general 
one could distinguish more complicated ordering structures too. They occur 
at the strong couplings, to be addressed in Section 5. 



We are going to explore the phase diagram by measuring in particular the 
order parameters 

((^o) and ((^i) . (2.15) 

Based on the magnitudes of these expectation values we distinguish three 
phases, as we specify in Table [H Equivalent tools were applied before in 
investigations of the 2d model on a fuzzy sphere [S]. 

• In the disordered phase {(pe) ~ holds for all i. The angular mode 
decomposition does not detect any contribution that could indicate 
a spontaneous breaking of the rotational symmetry on the sphere. 

• The uniform ordered phase is characterised by ~ ((/^q) ^ , 
i.e. the zero mode contributes significantly, whereas higher modes are 
suppressed. This phase corresponds to the spontaneous magnetisa- 
tion in a ferromagnet. 

• In the non-uniform ordered phase a non-zero mode condenses, which 
leads to the relation 

(ii*ir)»(v^o)^o. 

In this case the rotational symmetry of the sphere is spontaneously 
broken. For the settings to be explored in Sections 3 and 4, this is 
manifest by a dominant contribution for i = 1 : (||<l>||^) «i {(pf) S> 0, 
i^Ui) ^ . 

In the case of strong coupling the non-uniform ordered phase is dom- 
inated by the condensation of higher modes, (ipi) ^ {ipi) ~ {ipo) ~ 
for some i > 1. 

The general order parameter for this phase reads — (/^q). 
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For a precise identification of tlie pliase transition lines, we also consider 
the susceptibility-type observable^ 



(2.16) 



which display peaks at the corresponding phase transitions. 

To further substantiate the measurement of the phase diagram we take 
thermodynamic quantities into consideration as well, in particular the in- 
ternal energy U and the specific heat C, 



A peak in C , and in one of the susceptibilities Xi ? indicates a (regularised) 
second order phase transition. 

Note that the non-uniform ordered phase is specific to the fuzzy sphere; 
it does not occur in a regularisation on a sharp sphere, or — generally 
speaking — on commutative spaces. In the fiat non-commutative space, 
such a phase was predicted in Ref. [H] for the Xcj)"^ model in 3 and 4 
dimensions, as a consequence of the notorious mixing of ultraviolet and 
infrared singularities (UV/IR mixing). For d = 4 arguments involving the 
renormalisation group |T7] and an effective action [18] were added. In = 3 
this behaviour could in fact be demonstrated numerically [19] by means 
of lattice simulation results, which were extrapolated to a simultaneous 
UV and IR limit while keeping the non-commutativity constant ("double 
scaling limit"). We will discuss the relation between that result and the 
system studied here in Subsection 5.2. 

These three phases, including the phase of non-uniform order, were 
also observed numerically on the fuzzy sphere without time direction f8l|, 
in agreement with theoretical considerations [20^ i21j . Similarly this exotic 
phase was found in lattice studies of the 2d non-commutative plane j2^ITU] . 
However, in that double scaling limit has not been worked out so 

far, hence the existence of this phase in the continuous plane is an open 
questionJ§ 

Here we reconsider the 3d model. However, the spatial part is not 
accommodated on a non-commutative plane but on a fuzzy sphere, as we 
pointed out before. In addition our main interest refers to the extrapolation 
to a commutative limit — in contrast to double scaling limit addressed 
in Ref. [19] — in view of the possibility of using the fuzzy sphere as a 
regularisation scheme for ordinary (i.e. commutative) field theory. 

^In the following we will refer to them simply as "susceptibilities" . 

^Due to the non-locality it is not ruled out by the Mermin- Wagner Theorem. Still 
Ref. 16J does not expect this phase (for a charged scalar field) in rf = 2, based on an 
extension of this Theorem to a related effective action with an unusual kinetic term. 
For a neutral scalar field Ref. [23] arrives at the opposite conclusion. 



U = {S) , C= {S^) - (5)2 . 



(2.17) 
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We performed all our simulations at 



Nt = N , (2.18) 

so that the system has the same number of degrees of freedom in the 
temporal and in the spatial directions. 



3 Determination of the phase diagram 

To explore the phase diagram we fixed some value of A and varied 
searching for a phase transition. Decreasing is analogous to lowering 
the temperature in statistical mechanics. In all settings we could identify 
a critical value < 0. For > we are in the disordered phase 
(corresponding to high temperature), whereas < gives rise to the 
dominance of some ordering, and therefore spontaneous symmetry break- 
ing. For small values of A this order is uniform (like the spontaneous 
magnetisation of a ferromagnet), but for larger A it becomes non-uniform 
(some kind of staggered order), cf. Section 2. 

Let us describe the determination of m^. As a first example. Figure 
[1] shows results at = 16, i? = 4, A = 0.44. The specific heat takes its 
maximum at ^ —0.32 , and the susceptibilities confirm this value. The 
peak for xo further specifies that we enter the uniform ordered phase for 

< ml- Unlike C, the xe are sensitive to the type of order below m^. 



Xo I — ' — I 

i 



-0.4 -0.36 -0.32 -0.28 -0.4 -0.36 -0.32 -0.28 

2 2 

m m 

Figure 1: The specific heat C (on the left) and the susceptibilities xo ond 
Xi (on the right, defined in eq. Ii2.16\) ) for N = 16, R = 4 and A = 0.44. 
The location of the maximum of C coincides with the peak in xo? which 
provides a consistent result for the critical value m?^ ~ —0.32 . 

This agreement between the two criteria gives a reliable determination 
of ml- Figure [2] shows this consistency in the case = 12, i? = 8 for a 
variety of A values. It also gives an overview of the phase diagram: as A 
rises, ml moves to more negative values. This relation is linear to a good 
approximation, as we are going to discuss in Section 4. 
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Figure 2: The order/ disorder transition line in the phase diagram for N = 
12, R = 8. The transition is identified consistently from two criteria. This 
figure shows the regime of weak coupling, where m? < ml implies a uniform 
order. 

For fixed values of N and R, tlie (A,m^)-plane contains a triple point, 
which we denote as {XT^rn^). It separates the regimes of weak coupling, 
A < At, and of moderate or strong coupling, A > A^ or A ^ Ay. 

A typical example for the behaviour of the order parameters at weak 
coupling is shown in Figure [31 For sufficiently negative m? the order pa- 
rameter ((y9o) rises drastically. The peaks in xo and in the specific heat 
allow for a more accurate evaluation of ml = —0.12(2). Further insight 
into this phase transition is gained by splitting the internal energy U (in 
eq. 02.171) ) into contributions due to the different terms in the action (12.61) . 



U2 = f (TrE,(**-C^$0]) 
f/3 = ^(TrE,<f?]) 
^4 = 4^(TrE,$^]) 



spatial kinetic contribution 
temporal kinetic contribution 
contribution due to the mass term 
contribution due to the self-interaction. 



They fulfil the identities 

Ui + U2 + Us + U^ = U , 2(f/i + f/2 + f/3 + 2f/4) = 1 , (3.1) 

(the latter is obtained from a variational argument). In our setting the 
total internal energy U is identical to the entropy. The last plot in Figure 
E] illustrates that it deviates from a constant as is decreased below m^. 
Moreover we see that U3 and U4 drift away from zero at this transition: 
the quadratic term becomes negative and the quartic term positive, which 
is just the situation that triggers spontaneous symmetry breaking with 
uniform ground states. 
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Figure 3: Simulation results for the determination ofrn^ at N = 12, R = A, 
A = 0.17. On top we show the order parameters according to specification 
Ii2.15\) (on the left) and the corresponding susceptibilities (on the right). We 
observe a transition between disorder and uniform order at ml = —0.12(2). 
Below: the specific heat (on the left) and different contributions to the in- 
ternal energy U = J2t=i T^'^ right). In the phase of uniform order, 
the potential terms deviate significantly from zero, so that the mass term 
f/3 (the quartic term U4) becomes negative (positive). The kinetic contribu- 
tions Ui (spatial) and U2 (temporal) remain almost constant — here each 
dimension contributes approximately the same amount. 



Let us proceed to the regime of moderate coupling. Figure H] shows the 
order parameters for an example in that regime, along with the suscep- 
tibilities. The data reveal a transition between disorder and non-uniform 
order at ml = —0.37(2). To provide an overview, we sketch in FigureOthe 
complete phase diagram obtained at = 16, i? = 8, as a further example. 



4 The scaling of the phase transitions 

In this section we analyse the scaling of the observed phase transition lines 
— and in particular their intersection in the triple point — with respect 
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Figure 4: T/ie order parameters and susceptibilities for an example at mod- 
erate coupling, N = 12, R = 8, X = 1.25. At ~ —0.37 we observe a 
transition between disorder and non-uniform order. 





-0.1 
-0.2 
-0.3 
-0.4 
-0.5 
-0.6 
-0.7 



transition at weal< coupling 
transition at moderate coupiing 
triple point 



O/Vy 



X" 



0.5 



1.5 



Figure 5: Overview of the phase diagram at N = 16 and R = 8. 



to and R. We proceed by considering separately the boundaries of the 
disordered phase with the two ordered phases. 



4.1 Transition between the disordered and the uni- 
form ordered phase 

We first return to the fixed size = 12, i? = 8. Figure [2] shows the 
measured disorder /uniform order transition hne, 

ml = -0.39(2)A . (4.1) 

Probing also other sizes — by varying and R separately — we observed 
linear relations again, so we were guided to the ansatz 

ml = f{N, R) X , f{N, R) oc N^^'R^'' . (4.2) 
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In the last expression we anticipate that the function f{N,R) can be pa- 
rameterised successfully in a monomial form. To illustrate this, we first 
fix again R = 8 but vary in the range 8 ... 33. Figures [6] and [7] (on the 
left) show that we obtain excellent fits with the exponent 6n — 0.64 at 
i? = 4 ... 16. By including similar plots for R = 2 and 32 we extract 

Sn = 0.64(3) . (4.3) 

Subsequently we explore the dependence of f{N, R) on R, and we find 
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Figure 6: The phase transition disorder/uniform order at N = 8 .. .33 for 
R = 4 (on the left), and for R = 8 (on the right) . We observe consistency 
with the exponent Sp^ ~ 0.64 given in eq. ( [^.5] ). 

again agreement with the monomial ansatz (14.21) . Moreover, it turns out 
that the function / essentially just depends on the ratio N/R. The expo- 
nent 6r and the coefficient are determined from the fit in Figure [7] (on the 
right), and we arrive at 

^0.64(3) 

= -0-31(1)^^^ A , (A < At) . (4.4) 

We add that the quahty of the corresponding fits is very good as long as 
3/8 < N/ R<8. If the ratio N/ R leaves this interval on either side, the data 
request a more general function f{N, R) — beyond the monomial form — 
though the ansatz for TTT/^ clS db linear function of A remains applicable over 
a wider range. 

4.2 Transition between the disordered and the non- 
uniform ordered phase 

Figure [8] (on the left) shows the extension of Figure [2] to much stronger 
couplings A. The critical values now mark the transition between dis- 
order and non-uniform order. Over this range of A the phase transition 
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Figure 7: On the left: the transitions disorder/uniform order at R = 16 
for a variety of N values. On the right: the fit for the function fi{R) := 
f{N, R)N~'^-^'^, which yields eq. fi4-4\ ) for the phase transitions at weak cou- 
pling. 



line is curved. However, at this stage we concentrate on the vicinity of the 
triple pointjll which is around {Xt^tti^) ~ (0.60, —0.23) in this diagram. 
Up to A ~ 3 a linear fit for is again very precise, but it now requires an 
additive constant, ml = -0.11(2) - 0.22(1)A . 
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Figure 8: The transition disorder /non-uniform order at R = 8 and N = 12 
(on the left), and = 8 ... 33 (on the right). 

We extend also this consideration to = 8 . . . 33, see Figure [8] on the 
right. A global fit implies the generalised form 

ml = -0.125(25) - 0.0178(10)ArA . (4.5) 

In analogy to Subsection 4.1 we repeated this study for i? = 2, 4, 16 
and 32; the case i? = 16 is shown as an example in Figure (on the left). 
Since the choice of A^A as the parameter on the works well in all 

''Large A values will be addressed in Subsection 5.3. 
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these cases, we are led to the ansatz 

ml = -g{R) - h{R)NX , (A > At) . (4.6) 

Figure [9] (on the right) shows our fits for the functions g{R) and h{R), 
which imply 

g{R) = 6.0(15)i?^i-^2(io) ^ ^ 0.082(7)i?-°-^2{4) _ ^4 7^, 
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Figure 9: On the left: the transition disorder/non-uniform order in the 
phase diagram at R = 16. On the right: the fits for the functions g{R) and 
h{R), which lead to the form ( [.^.71 ). 



4.3 The triple point 

We are now in a position to evaluate the triple point by the intersection of 
the phase transition lines at low and moderate coupling, given in eq. (14.41) . 
and in eqs. (14. 6p . (14. 7p . The intersection point is located at 

(l9(5)i?-i-20{ii)iV-o-64(9)^ -5.9(16)i?-o-56(")iVO-o(i)) 
(AT,my)~ /^o.o8{4) _ o.26(2)Aro-36(9) ' ^^'^^ 

We assume Gaussian error propagation, and include a possible error on the 
power of in the ansatz (14. 6p . 

This equation parameterises the triple point in the range for and 
R that we simulated. However, this formula cannot be extrapolated to 
large A^, due to the form of the denominator]^ A direct investigation of the 
triple point — and therefore of the phase diagram — in the thermodynamic 
limit would require simulations at very large A^. We will come back to this 
issue in Section 6, where we conjecture the properties of this extrapolation 
indirectly. 



^Unless one fixes i? (x iV, such extrapolations also takes us beyond the validity 
interval for the ratio N/R that we specified in Subsection 4.1. 
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5 Confrontation with related models 



Having elaborated the features of the phase diagram, we now compare it 
to related models, which were also simulated in the recent years, and to 
reduced matrix models, which have been studied analytically. 

5.1 The 2d model on a fuzzy sphere 

At finite A^, the phase diagram for the 2d model on a fuzzy sphere is 
qualitatively the same as we found in Section 3. An explicit convergence of 
our model to the 2d case could be expected in a setting which renders solely 
the temporal kinetic term negligible. We saw, however, in Figure [3] that 
the impact of this term remains significant and approximately constant at 
weak and moderate coupling. This property is generic in our study, hence 
there is no basis for expecting a reduction to the 2d model on a fuzzy 
sphere. At large A some reduction sets in, but it is of a different kind, see 
Subsection 5.3. 

5.2 Comparison with a non-commutative torus 

Ref. jl9] presented a related numerical study of the 3d Xcj)'^ model: also 
there the Euclidean time was lattice discretised, while the spatial dimen- 
sions were treated by non-commutative coordinates. However, a lattice 
formulation and periodic boundary conditions were assumed for the spa- 
tial directions as well, and the non-commutativity tensor G^j, was constant. 
By means of Morita equivalence [23] the scalar field — first defined on a 
N X N lattice in the non-commutative plane — was mapped onto Hermi- 
tian N X N matrices, so the action ultimately simulated was similar to eq. 
(12.61) (and also the rule (12.181) was the same). The difference is the form of 
the spatial kinetic term: on the non-commutative plane it was constructed 
by an adjoint matrix operation which represents a shift by one lattice unit. 

The phase diagram at finite A^ was qualitatively equivalent to the form 
that we found here. On the torus the non-uniform ordered phase was de- 
noted as "striped phase", and the stripe formation for the two signs of 
0j(xi,X2) could indeed be visualised by mapping the matrices back to 2d 
lattice configurations (at a fixed time site). 

Let us proceed to a quantitative comparison. On the torus the boundary 
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of the disordered phase was approximately identified as [IS] 



disorder / uniform 
disorder/ non— uniform 



triple point 



ml ~ -0.80A 
ml ~ -0.48A - 



80 

iV2 



(250, -200) 



(5.1) 
(5.2) 
(5.3) 



The transition disorder /non-uniform order was only explored up to mod- 
erate coupling, A = 0(1). 

Our corresponding results are given in eqs. (14. 4p and (14.61) to (14.81) . The 
transition line (15.11) can be matched if we set R ~ 0.23A^, but the triple 
point still differs strongly. However, in the weak coupling region a link 
between the two models can hardly be expected, due to the significance of 
the spatial kinetic term. 

We proceed to moderate coupling and look at the requirements for 
agreement with eq. (15.21) . The condition for the function g{N,R) reads 
R ~ 0.26iV^ °^, which is close to the above requirement for eq. (15. ip . On 
the other hand, the condition due to the function h{N, R) deviates more 
from this pattern, in particular in view of the exponent {R ~ 0.28N^-^^). 

We anticipate at this point that the geometrical picture — to be de- 
scribed in Section 6 — suggests that the non-commutative plane emerges 
for 



Hence the value for [3 obtained at moderate A takes us even further away 
from the geometrical picture. This may appear somewhat surprising, but it 
is not paradoxical. For A/At = 0{1) the spatial kinetic term is significant, 
as the example in Figure [3] shows, so a difference in this term may well 
displace the phase transition lines. 

In addition we saw that our results in Section 4 suffer from finite 
artifacts, hence the considerations in this and the previous Subsection have 
to be interpreted cautiously. 

5.3 Reduction to a matrix chain 

At strong coupling A the impact of the kinetic terms is reduced. In the 
extreme case where they are fully negligible, one could imagine a transition 
to a 1-matrix model of a single Hermitian random matrix $i with the 
potential 47ri?2-Tr [^<l>2 + |$4] . That model has been studied analytically 
at > [25] and at < with the result [26] 



R(xN^ , (3=1/2 . 



(5.4) 




(5.5) 
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Consistent numerical data have been reported for the 2d model on a 
fuzzy sphere [S] and on a non-commutative plane [27j. In our simulations 
we explored values up to A = 0(10^ ... 10^), but we could not observe the 
feature of eq. (15.51) [llj. 

We did find, however, consistency with a partial reduction, which only 
neglects the spatial kinetic term; this is opposite to the scenario commented 
on in Subsection 5.lJ^ In our formulation (12. 6p this means that the double- 
commutators [Ci, [Ci, are negligible. This behaviour is exactly con- 
firmed by Figure [TOl Thus we obtain a matrix chain (analogous to a spin 
chain), consisting of Hermitian matrices with a quartic potential, which 
are linked by a discrete second derivative. In the large N limit Ref. f2E\ 
derived the critical line 



This formula corresponds to the one-cut vs. two-cut transition in the eigen- 
value distribution of random matrices. This transition plays an essential 
role; it must show up at least for sufficiently large A. In addition it is of 
importance down to the triple point, where it merges with the transition to 
the uniform order. But as the coupling is lowered, the fuzzy kinetic term 
distorts the form (15. 6p for this transition line. It eventually terminates at 
the triple point — a feature not present in the matrix chain. 

Figure [TT] shows two examples where the behaviour of eq. (15. 6p is 
matched accurately, up to a modification in the coefficient. We take a 
closer look at the first case, = 8, i? = 16, where the precise agreement 
with the exponent in eq. (15. 6p is amazing because of the relatively small 
matrices. We further fix A = 625 and we show in Figure [12] the specific 
heat and the order parameters. The former confirms the critical parameter 
ml = —9.0(5) which appears in Figure [TTl The order parameters demon- 
strate that we deal with a transition between disorder and non-uniform or- 
der, where the latter corresponds now to the condensation of higher modes, 
i > 1 (unlike the examples in Subsection 4.2). 

This type of model (matrix quantum mechanics) has applications in var- 
ious branches of physics. For instance QCD with Nf light quarks fiavours 
in a small box but an elongated time direction (the so-called "5-regime") 
can be treated effectively by quantum mechanics of SU{Nf) matrices [28] . 
The very same technique can also be applied in solid state physics |29j . 

Our case of Hermitian matrices attracted attention in string theory 
since the early nineties [SUl [SI] • In that framework it represents the c = 1- 
model, which describes random surfaces moving in one dimensions. At 

^One should expect the model on a non-commutative torus to agree in this regime, 
but it has only been simulated up to moderate coupling |19j . as we mentioned in Sub- 
section 5.2. 




(5.6) 
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Figure 10: The different contributions to the internal energy at N = 8, 
R = IQ and A = 625. The temporal kinetic contribution U2 now clearly 
dominates over the spatial kinetic contribution Ui, in contrast to the weak 
coupling behaviour shown in the last plot of Figure O This is fully consis- 
tent with the observed reduction in the spatial directions only, which leads 
to the matrix chain behaviour described in Ref. IM 
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Figure 11: The phase transition at strong coupling X for N = 8, R = 16 
(on the left), and for N = 23, R = 16 (on the right). In both cases we 
observe a broad window of agreement with the matrix chain formula / I5. 6\) . 



a finite time periodicity T there is a vortices-driven Kosterlitz-Thouless 
phase transition to the c = 0-model [30l [32] . Possible hnks to QCD2, to 
2d black holes and to topological field theory were studied intensively [33] . 
This model continues to attract interest, see Refs. [3l] for recent examples. 
An overview of the different reduction scenarios is added in Table [2l 
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Figure 12: The specific heat (on the left) and the order parameters (on 
the right) at N = 8, R = 16 and A = 625. We confirm the critical value 
= —9.0(5) at the transition between disorder and a non-uniform order 
(cf. Figure [71] on the left). The plot on the right demonstrates that it 
corresponds to a higher mode condensation, i.e. to a non-uniform order 
with a complicated pattern. 



spatial 


temporal 


setting 


status 


kin. term 


kin. term 










2d model on 


numerical studies [8] 


large 


small 


a fuzzy sphere 


not attained here 






matrix chain 


analytical studies [251 126j 


small 


large 


or c=l-model 


reproduced here at large A 






total reduction 


analytical study [26\ 


small 


small 


1-matrix model 


not attained here 



Table 2: An overview of the conceivable reductions of our system — due to 
the negligibility of kinetic terms — and their status in the literature and in 
this work. 

6 A conjecture about the large N limit 

In this section we discuss the extrapolation of the phase diagram to large 
(which represents the thermodynamic limit) and to large R (the transition 
of the spatial part to a plane). In addition our rule (12.181) connects the 
large N limit with the thermodynamic limit in the temporal direction. 

We first consider the geometry of the sphere under these limits. As 
we remove the cutoff A^, the coordinates (12. 3p describe different 2d spaces, 
depending on the simultaneous treatment of the radius R: 

• The limit N —>■ oo at R = const, leads to a sharp sphere. 
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• If the radius grows slowly, R oc A^^, 0</3<l/2, we end up with a 
sharp plane in the large N limit. 

• If we take instead oc i?^ — oo we obtain a non-commutative 
plane with a constant non-commutativity tensor Q^i, = i 6 e^p (where 
/X, 1/ e {1,2}). This can be seen for instance in the plane which 
emerges around the point (0, 0, 1), 

[^i,^2]^i^2£L^ ^ e = '^-^. (6.1) 

ViV2 - 1 R N ^ ' 

The scaling behaviour of a field theory on this space, and in particular 
its phase diagram, still has to be investigated. For the parameter range 
simulated, this was carried out in Section 4. We now address the issue 
of the large N extrapolation, which was postponed in Section 4. For the 
disorder /non- uniform order transition, we found agreement with the eqs. 
(14.61) and (14. 7p at moderate A, and with eq. (15.61) at strong A (up to a 
modest modification of the coefficient). In either case A only occurs in 
a product with a positive power of (1 resp. 4/3). This is consistent 
with the suppression of field fluctuations at large or large A. Based on 
this property, we conjecture that exploring the large A^ behaviour could be 
equivalent to the case of large A at the values of A^ in our study. 

This means that we now refer to eq. ( 15. 6p to determine the triple point. 
The behaviour 



m 



2 



-c(Ar/i?)4/3A2/3 (c = const.) (6.2) 

is compatible with our data. For instance at i? = 16 we obtained m^A~^/^ ~ 
-0.121 at A^ = 8, and -0.31 at A^ = 16 (see Figures [TT] and [T3l) . which 
matches very well the behaviour ml oc N'^/'^ . Also the dependence of the 
radius follows eq. (16.21) : for example the results at A^ = 8 and = 8 vs. 
-R = 16 (in Figures [TT] and [T3|) agree very well with the relation ml oc R~^^^. 
Considering many flts of that kind [TT] we found that both exponents 
and the coefficients c in eq. (16. 2p fluctuate within about 10% around the 
theortical values of eq. (15.61) . which were derived for the matrix chain at 
large A^ [26]. 

The intersection of this curve with the disorder/uniform order transition 
line yields 

jY2.1(2) 

OC . (6.3) 



In the framework of this conjecture, we obtain the following scenario^'^^' 
for the limit N oc: 



-"^"Based on the errors of the exponents in eq. (|6.3[) . the distinction between the scenar- 
ios should actuaUy refer to /? = 1.0(1) (where /3 is defined in eq. (|5.4[) ). but for simphcity 
we just refer to (3 — 1. 
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Figure 13: The phase transition at strong coupling A for N = 16, R = 16 
(on the left), and for N = 8, R = 8 (on the right). As in Figure [771 we 
find agreement with the matrix chain prediction ml oc A^/^. In addition, 
comparison to the result for N = 8, R = 16 in Figure [771 supports the full 
proportionality relation ( fg.^j . 



• Limits with oo, N/R — > oo remove the phase of non-uniform 
order. Then the ordered regime only consists of the uniform phase, 
which is separated by an Ising transition from the disordered phase. 
This hmit corresponds to a commutative modeL It includes in partic- 
ular the case of a sharp sphere {R = const.). In contrast to the tree- 
level expectation, this class of limits also captures the case N oc R^, 
which geometrically leads to a non-commutative plane. 

• If the limits N, R ^ oo are taken such that the ratio N/R remains 
finite, the triple point stabilises and the phase diagram keeps qualita- 
tively the form that we observed at finite A^; all three phases persist. 

• Limits with oo, N/R — > remove the phase of uniform order. 
Now the ordered regime consists solely of the non-uniform phase. 
This scenario is obtained for a rapidly expanding sphere. Here the 
non-commutativity dominates the thermodynamic limit. 

The transition lines that we identified at small and at moderate coupling 
in Section 4 do not scale simultaneously for for any fixed choice of the axes 
of the phase diagram, as the forms (14. 4p and (14.6!) show. On the other 
hand, the large N behaviour (as conjectured in this section) overcomes 
this problem: if the triple point moves to or oo, only one transition line 
survives. In the case R cc N, which stabilises a finite triple point, the axes 
(A,m^) apply to both transition lines, without the necessity of rescaling. 
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7 Conclusions 



We presented a numerical study of the phase diagram in the 3d Xcp^ model, 
where the spatial part is regularised on a fuzzy sphere, while the Euclidean 
time is lattice discretised. On the regularised level, we identified three 
phases based on the order parameters ipe in eq. (12 .14^ . the corresponding 
susceptibilities and the specific heat. At fixed A there is a critical parameter 

< : for > (m^ < m^) the system is disordered (ordered). The 
transition to < leads to a uniform order at small A, and to a non- 
uniform order at moderate or large a1"1 

The boundary between these two scenarios corresponds to the triple 
point, which we denoted as {XT,'m'^). The transition disorder /uniform 
order (at A < At) is analogous to a spontaneous magnetisation and its 
critical line is parameterised by eq. (14. 4p . 

The non-uniform ordered phase emerges as a consequence of the non- 
locality in the fuzzy sphere regularisation. That phase does not occur in 
a pure lattice regularisation. It corresponds to a spontaneously broken 
rotation symmetry. At moderate coupling strength, A > At, the critical 
line is described by eqs. (14. 6 p and (14.71) . From its intersection with the 
curve (14. 4p we infer the location of the triple point given in eq. (14. 8p . This 
formula captures a large amount of data that we collected, but it cannot 
be extrapolated to the limit N oo. 

Next we discussed the relation between the system studied here and 
some models investigated in the literature, which have qualitatively similar 
phase diagrams at finite N. We found, however, significant differences from 
the 2d A0^ model an a fuzzy sphere [S] . As for the 3d A0^ model on a non- 
commutative torus [in] the vicinity of the triple point can be matched 
rather roughly for a suitable relation between and R. 

At large A, we observed consistency with the reduction to a matrix 
chain model, which had been solved in the large A^ limit [26] . In particular 
our data match precisely the predicted relation oc {N'^X/R'^Y^^. That 
model is known as the c = 1-model in string theory [SUl ED E21 E31 Elj, 
which we have therefore captured non-perturbatively. 

We then conjectured that large A^ values may be equivalent to large A 
couplings at moderate A^, since A^ and A tend to appear only as products 
in the formulae for the phase transition lines. This conjecture leads to 
different limits depending on the exponent f3 in the relation R oc A^^. For 
/? < 1 we obtain a commutative limit, with an Ising-type transition between 

^^We always refer to a region where rri^ is kept of the same magnitude as m^; driving 
it to <C causes simulation problems with the thermalisation and decorrelation, 
hence we could not explore that region reliably. Similar technical problems obstructed 
a direct observation of the transition between the two ordered phases. 
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disorder and uniform order. On the other hand, a rapidly expanding sphere 
{(3 > 1) leads to a dominance of the non-uniform phase. That feature is 
a characteristic for a non-commutative theory, where UV/IR mixing gives 
rise to an ordering due to the condensation of a non-zero mode |T6l [T71 [T8|, 

m\. 

It remains an open question to verify the conjectured limits by direct 
inspection of the triple point at very large system sizes N. For instance, 
at R = const, one should verify the properties of the fuzzy sphere at larger 
A^. In the light of gravity- induced non-commutativity on the Planck scale 
[35] , a relation to quantum effects in a black hole might be conceivable (see 
e.g. Refs. [221 ES] for this line of thought). 

According to our conjecture, the distinction between the scenarios of a 
commutative and a non-commutative limit is not located at the point where 
it is expected on purely geometrical grounds (/? = 1/2). Furthermore, it 
does not coincide with the picture suggested by a perturbative calculation 
to two loops [20]. In that picture the commutative continuum limit could 
not be retrieved at all because of the UV/IR mixing [2D]@ 

A general lesson is that non-locality — once it is introduced — can 
cause surprises in the extrapolations on the non-perturbative level. This 
observation may serve as a warning also for the use of non-local actions in 
lattice simulations, such as rooted staggered fermions (cf. footnote Hj) or 
overlap fermions [STj at strong gauge coupling0 On the other hand, such 
surprises may provide insight into other universality classes of interest, as 
we have seen. However, using the fuzzy sphere as a regularisation scheme 
in quantum field theory is not straightforward and its application requires 
a careful investigation of the phase diagram. 

A Technical aspects of the simulation 

Our simulations were based on the Metropolis algorithm. In each step we 
updated only one pair of conjugate matrix elements, {^t)ij = {^t)*ji ■ A 
detailed comparison in technically similar simulations revealed that this is 
far more efficient than updating complete matrices $f [T9] . It proved useful 
to propose independent changes of the real and the imaginary part of {^t)ij 
with absolute values below N ^Jm? /\ (and flat probability distribution); 

^^Thc model studied perturbatively in Refs. [20] coincides with the one considered 
here, up to the use of a continuous Euchdean time. In the framework of the c — 1- 
model (at large A) it is a mystery that this difference seems to imply a T ~ 1 /T duality, 
which is absent in the matrix chain [30 . 

^•^Despite the inverse square root in the overlap operator, it is local at weak gauge 
coupling resp. on fine lattices [55] ■ The range of locality — and therefore of a safe 
definition of chiral fermions — can still be enlarged by a non-standard kernel [32] ■ 
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this led to acceptance rates typically around 1/2. One sweep applies this 
step successively to all independent elements in a configuration $. 

We were often confronted with several local minima of the action. In 
many cases these minima were too pronounced for a tunnelling to occur 
even in histories involving O(IO^) sweeps. This property obviously ob- 
structs the direct measurement of observables. As a first remedy we per- 
formed a number of runs with independent hot starts. At the end we 
summed up the statistics collected in each run after thermalisation (which 
we are going to comment on below). The histories in all runs had the same 
length of O(IO^) sweeps. This improves the situation, but the number of 
runs was still too small (typically 0(10)) to sample the vicinities of the 
different minima reliably. 

Therefore we extended the algorithm as follows. We stored the end 
configuration of each run as $end- The subsequent run takes a new hot 
start, but after thermalisation the current configuration may be replaced 
by $end through a Metropolis accept/reject step, before the history contin- 
ues. We denote this method as adaptive Metropolis algorithm. In fact it 
improves the statistically correct inclusion of the vicinities around various 
minima, and it leads to stable and sensible measurements, as the examples 
in Figure [T3] illustrate. The impact of higher local minima tends to be 
overestimated by fully independent runs. The additional Metropolis step 
helps to overcome this artifact; in particular the global minimum now re- 
ceives the suitable weight. Figure [13 (on the left) shows an example for 
this effect 

In principle the variety of metastable vacua can be regarded as a severe 
thermalisation problem. However, since it is taken care of by the adaptive 
Metropolis step, we reduce our notion of thermalisation to the Monte Carlo 
time in each run until the observables stabilise over a long period at the 
value that corresponds to the chosen minimum. In this respect, about 2000 
sweeps were sufficient to thermalise quantities like the action and our order 
parameters (introduced in Section 2). 

While this thermalisation is harmless, a technical problem could oc- 
cur due to a large number of local minima, in particular at A ^ At- 
For a simplified consideration we assume the kinetic terms to be negli- 
gible (although we saw in Subsection 5.3 that this complete reduction is 
not really achieved). Then the matrices $t are independent and — in a 
minimum of the potential — each one can be transformed to a diagonal 

^"^In his study of the 2d model on a fuzzy sphere, M. Panero apphed successfully an 
"overrelaxation" technique, which is described in his works quoted in Ref. [8]. However, 
this technique is unlikely to be applicable in our case, due to the presence of the temporal 
kinetic term. 
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Figure 14: An illustration of the progress due to the adaptive Metropolis 
algorithm. The parameters in these examples are N = 12, R = 4, X = 1.83. 
On the left we show the specific heat C ; for adaptive Metropolis only the 
peak is clearly observed at the critical value m^. The plot on the right is 
a histogram for Tr(<i>) at m? = —1.7, i.e. in the phase of uniform order. 
The two peaks approximate the expected symmetric form well with adaptive 
Metropolis, but not with the same number of independent runs. 

form with elements ±^y\m'^\/(NX). With all sign combinations the term 
Coo = Tr(<l>)/A^^ can take N'^ + 1 values, which is a considerable num- 
ber for the system sizes that we studied. However, for independent random 
signs the values near zero dominate, whereas the probabilities for minima 
with large ipo = |coo| are suppressed. For instance Figure [15] (on the right) 
shows a histogram for Cqo at = 12 ; only 5 peaks (corresponding to the 
5 dominant minima) are visible. 

1 1 1 1 1 1 1 I — I — I — I — I — I — I — I — I — I — 

independent runs 




Figure 15: Histograms for two quantities, which are sensitive to the variety 
of action minima. On the left: the internal energy U measured at N = 12, 
R = A, \ = 1.83, m2 = -1.7. On the right: Cqo (see eqs. [2ll\) . [2JE) ) 
measured at N = 12, R = 16, A = 0.75, = —0.3. Coo is shown in 
units of 2^ 'n:\vrfi\l (N^I'^X) ~ 0.0045 so that the peaks are located at integer 
values, in agreement with our prediction. 

The statistical errors were evaluated independently with the binning 
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and the jackknife method on one hand (we probed various bin sizes), and 
with the Madras-Sokal method [ID] on the other hand. The latter amphfies 
the standard error by a factor which takes the autocorrelation into account. 
We generally display the largest (and therefore safest) error bar obtained 
by these methods. 
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